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ABSTRACT 

Galaxy clusters are not iso-thermal, and the radial temperature dependence will affect 
the cluster parameters derived through the observation of the Sunyaev-Zeldovich (SZ) 
effect. We show that the derived peculiar velocity will be systematically shifted by 
10 — 20%. For future all-sky surveys one cannot rely on the observationally expensive 
X-ray observations to remove this systematic error, but one should instead reach for 
sufficient angular resolution to perform a deprojection in the SZ spectra. The Compton 
weighted electron temperature is accurately derived through SZ observations. 
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1 INTRODUCTION 



Galaxy clusters have been known and studied for many 
years, and the radial dependence of cluster temperatures is 
becoming a testing ground for models of structure formation 
and for our understanding of gas dynamics. Galaxy clusters 
are not iso-thermal. Instead, the emerging temperature pro- 
file is one where the temperature is approximately flat or 
increases from the centre to some characteristic radius, and 
then decreases again for larger radii. 

The central temperature decrement has been much dis- 
cussed and the possibility of cooling flows has been explained 
in excellent reviews (Fabian 1994; Donahue & Voit 2003). 
Many clusters are well fit with a power law T ~ r T in the 
very central region (Voigt & Fabian 2003), with a slope, 
t, between 0.15 and 0.45. Numerical simulations are only 
now beginning to see this central decrement (see Motl et 
al. (2003) for references). 

The outer temperature decrease is well established 
both observationally (Markevitch et al. 1998; De Grandi & 
Molendi 2002; Kaastra et al. 2003; Arnaud, Pratt & Pointe- 
couteau 2003) and numerically (see Lin et al. (2003) and 
references therein). The simulations have even started to 
provide estimates of the outer temperature decrement which 
are well fit with power laws in fair agreement with observa- 
tions. 

Such non-trivial temperature profiles also affect derived 
cosmological parameters like the Hubble parameter (Battis- 
telli et al. 2003; Lin, Norman & Bryan 2003). 

On the other hand, observations of the Sunyaev- 
Zeldovich (SZ) effect are becoming increasingly accu- 
rate (Laroque et al. 2002; De Petris et al. 2002; Battistelli et 
al. 2003), however, most analyses of SZ observational data 
are made under the simplifying assumption of iso-thermality. 



Several groups have considered the ability of future SZ ob- 
servations to extract cluster parameters (Knox, Holder & 
Church 2003; Aghanim, Hansen & Lagache 2004), but al- 
ways under the assumption of iso-thermality. We therefore 
set out to study the effect on the SZ derived parameters of 
non-trivial cluster temperature profiles. 

In section 2 we show that the relevant quantities for 
SZ observations are Compton-averaged, in particular we 
emphasize the difference between the temperatures derived 
through X-ray and SZ observation for clusters which are not 
iso-thermal. In section 3 we use observed and simulated clus- 
ter profiles to show that the systematic shift in the derived 
peculiar velocity, v p , is of the order 10 — 20%. Finally we 
comment on the possibility of deprojection of SZ observed 
spectra in section 4. 



2 THE SZ EFFECT 

The SZ effect is dominated by the Compton parameter 
kT e 



y = o T 



■dl. 



(1) 



where gt is the Thomson cross section, n e ,T e and m e are 
number density, temperature and mass of the electrons, and 
the integral is along the line of sight. If we extract the ra- 
dial dependence of the parameters (Birkinshaw, Hughes & 
Arnaud 1991) 



n e (r) 
Te(r) 



n° fn e (r) , 
T°f Te (r) : 



then y can be written as y = k v J fn,,fT c dl, where the con- 
: (JTn° e kTe I (m e c 2 ) . It will be convenient to in- 
fn e /t c , in which case we have 
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Ky J ^ v 



dl. 



Similarly, the optical depth is 



(2) 



(3) 



Al 



with k t = cttu®, and one sees that for an iso-thermal cluster, 
where fr e (r) = 1, one has y ~ rT e °. 

For any cluster observation one receives information 
from an integral along the line of sight. Since the radial 
structures of clusters are non-trivial this implies that one 
may have to define averaged quantities. For X-ray observa- 
tions the averaging procedure is often simplified to emission- 
weighted quantities, scaled by the emissivity which for large 
temperatures roughly scales (Sarazin 1986) as e ~ r?T}J 2 



(Oh 



J Oedl 

J7dT 



(4) 



For the SZ effect the relevant scale is not the emissivity, 
but instead the local Compton parameter, n e (r)T e (r), and 
therefore the average of O is given by 



(0)y = 



S°fv 



Jfydl ' 

for instance, the average SZ temperature for a cluster is 

J T e fy dl Ky J Tefy dl 
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(6) 



We will hereafter only discuss Compton-averaged quantities, 
and hence omit the index on averages. 

The SZ effect for an iso-thermal cluster is composed of 
several parts 

^ = y (g(x) + T e S(x,T e )) - j3 p r a(x) , (7) 

where /3 P = v p /c is the bulk motion of the cluster referred 
to as the peculiar velocity. The spectral forms due to the 
thermal up-scattering, g(x), and due to the overall cluster 
motion, a(x), are well-known (Sunyaev & Zel'dovich 1972; 
Sunyaev & Zel'dovich 1980), and the spectral form of the rel- 
ativistic corrections to the thermal up-scattering, 5(x,T e ), 
is easily calculated (Wright 1979; Rephaeli 1995; Dolgov et 
al. 2001; Itoh & Nozawa 2003; Shimon & Rephaeli 2004). 
To a very good approximation 6(x,T e ) can be taken as 
independent of temperature (see e.g. Diego et al. (2002)), 
that is, for a given frequency, x, one can write 5(x,T e ) = 
bi(x) + bi(x)T R + 0(T 2 ), where 62 and higher order terms 
are subdominant. For the present discussion we can ig- 
nore a contribution to eq. (7) of the marginally detectable 
ultra-relativistic electrons (Ensslin & Kaiser 2000; Ensslin & 
Hansen 2004), since our findings will not depend on them. 
For excellent reviews on the SZ effect see (Birkinshaw 1999; 
Carlstrom, Holder & Reese 2002). 

For a cluster with a non-trivial temperature profile 
eq. (7) is written as an integral along the line of sight 



(8) 



~J~ = J [ Ky f y i 9 ^ + T "/T e 5(a;)) - /3pK T f„ e a(x)] dl , 

where the quantities f„ e (r) and f Te (r) are the local quanti- 
ties. Now, using the definition in eq. (5), this means 



y [g(x) + (T e )S(x) - (3 p {l/T e ) (m e c 2 /k) a(x)] . (9) 



For an iso-thermal cluster, where fr e = 1, one has l/(T e ) = 
(1/T e ), but for a thermally non-trivial cluster there are really 
4 independent variables, y,/3 p , (T e ) and (1/T e ). Clearly, for 
such a cluster the X-ray derived temperature, Tx , may differ 
significantly from (T e ) (as well as, 1/Tx =fc (1/T e )), and 
hence the X-ray derived temperature should only be used 
with great caution in the study of SZ observations. 

Most studies of SZ observations include up to 2 free 
parameters (Laroque et al. 2002; De Petris et al. 2002; Bat- 
tistelli et al. 2003) , for instance the Compton parameter and 
the peculiar velocity, and only very few have attempted stud- 
ies with 3 parameters (Hansen, Pastor & Semikoz 2002) with 
inclusion of the temperature. We will here consider the ef- 
fect of neglecting the non-trivial radial temperature profile, 
and show that it induces a systematic error on the peculiar 
velocity of the order 10%. 



3 CLUSTER PROFILES 

As a specific example we use a simple cluster model, where 
the electron density profile follows a /3-model 

n e (r)=n° e U + (£) J , (10) 

where r c is a characteristic radius. Such broken power-law 
density profiles are often used to fit observations (Cavaliere 
& Fusco-Femiano 1976; Lauer et al. 1995; Lewis, Buote & 
Stocke 2003), and are possibly even understood theoret- 
ically (Hansen & Stadel 2003; Hansen 2003; Ruszkowski, 
Bruggen & Begelman 2004), but may need slight modifica- 
tions for non- isothermal clusters (Ettori 2000) . The temper- 
ature profile may have the shape 



T e (r 



(11) 



where we use results from recent numerically simulated 
clusters (Lin, Norman & Bryan 2003) with a — 0.56 and 
P — 0.61. For cooling flow clusters the form is slightly more 
complicated, but an inclusion of the central decrement is 
straight forward. With these profiles one finds the average 
temperature from eq. (6) to be (T e ) = 7.08 keV when the 
central temperature is 10 keV, and (l/T e )(T e ) = 1.12 which 
indicates that one will find a peculiar velocity which is ap- 
proximately 12% systematically too large when using SZ ob- 
servations to define the cluster parameters. We will now test 
these findings numerically, by constructing and analysing 
the true SZ signal from this cluster model. 

With this cluster model, where we normalize the op- 
tical depth to r = 0.01, we construct the SZ signal along 
the line of sight through the cluster centre. As a specific 
example we take an observation with 4 observing frequen- 
cies at 90, 150, 220 and 270 GHz, and assume a sensitivity 
of \\iK for each channel. We choose v p = 1000 km/sec. For 
optimal observing frequencies (Aghanim et al. 2003; Holder 
2002) the statistical error-bars can be reduced slightly (rep- 
resented by A\ 2 = 1 on figure 1), but the systematic shift 
(about 12%) will remain the same. The systematic shift is 
independent of the assumed sensitivity, which here is taken 
slightly optimistic (ACT expects a sensitivity of 2/j.K). 
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Relative error 



Figure 1. The relative error on temperature, (T e ), and pecu- 
liar velocity, v p , through observation of a model SZ cluster with 
non-trivial temperature profile (eq. (11) with a = 0.56) and a 
/3-profile for the electron density (cq. 10 with (3 = 0.61). The SZ- 
derived temperature (solid, red line) is very close to the Comp- 
ton weighted temperature, but the peculiar velocity (dashed, blue 
line) is systematically shifted by 10 — 12%. 

One can now treat this constructed SZ signal (and 
corresponding error-bars) as a real observation, and de- 
rive the 3 cluster parameters, y,T e ,v v in the standard way. 
This is simply done by inputting the constructed SZ signal 
into the publicly available * SZ parameter extraction code 
sasz (Hansen 2004). The result is shown in figure 1, where 
we plot both peculiar velocity and temperature on the same 
figure. The solid (red) line is the relative error on the tem- 
perature, (T e — 7.08)/7.08, and the dashed (blue) line is the 
relative error on the peculiar velocity, (v p — 1000)/1000. As 
is clear, the average temperature is indeed very near the 
expected 7 keV, whereas the peculiar velocity is shifted by 
10 - 12%. 

The fact that the SZ determined temperature indeed 
turns out to be the "correct" Compton weighted tempera- 
ture implies that accurate prior temperature knowledge (e.g. 
from X-ray observations) would not change the systematic 
shift of the peculiar velocity. Having both accurate tempe- 
rature and electron density profiles from X-ray observations 
would allow one to calculate the systematic shift (ignoring 
the complications from an unknown clumping). However, 
peculiar velocities are most important when derived for a 
large sample of clusters (e.g. from an all-sky SZ survey), 
but detailed X-ray maps are observationally expensive and 
cannot be made for the large number of clusters expected in 
future all-sky SZ surveys. 

The minimal deviation from the derived central value 
of the temperature arises because S(x) in eq. (8) is not com- 
pletely independent of the temperature, but this minor effect 

* http: / /krone .physik.unizh. ch/~hansen/sz/ 

This SZ parameter extraction code uses the Monte Carlo like 
method of simulated annealing to derive the central SZ parame- 
ters and corresponding error-bars. 



is not important for our discussion. Furthermore, if one im- 
poses a "known" temperature of T e = Tx ~ 10 keV (from 
X-ray observations) then that would lead to a ~ 40% sys- 
tematic shift of v p , and even a "known" emission weighted 
temperature of Tx ~ 8 keV would give a ~ 20% systematic 
error on v p . To be very explicit, this means that if one knows 
the exact temperature and electron density profiles from X- 
ray observations, and then uses the emission weighted tem- 
perature in the analysis of SZ observations, then one would 
overestimate the peculiar velocity with 20% for this par- 
ticular cluster. If a given observation only puts bounds on 
the peculiar velocity, then the corresponding error-bars are 
overestimated by the same amount. This shows, as men- 
tioned earlier, that X-ray observed temperatures should be 
use with great caution in the study of SZ observations. If 
one does know the exact temperature and density profiles, 
then clearly one can calculate (1/T e ), with which one can 
extract the peculiar velocity correctly. 

We thus see that for this specific choice of temperature 
and electron density profiles the systematic shift in the de- 
rived peculiar velocity from SZ observations is 10 — 12%. 
The question now arises: which range in systematic shift is 
expected for real clusters? 

The density profiles of many clusters are well described 
by a /3-model with f3 ~ 0.67, however, numerical simulations 
show a rather large scatter in the numerical value for f3. 
Fixing the temperature profile to the one in eq. (11) with 
a = 0.56, and letting /3 vary from (3 = 0.75 to f) — 0.46 (Lin, 
Norman & Bryan 2003) leads to systematic shifts in the 
peculiar velocity of 8% to 20%. One could instead fix the 
density slope to /3 = 2/3, and use a polytropic equation of 
state 

with 7 in the range 1 to 5/3. Even though the polytropic 
shape does not provide an excellent fit to observations (De 
Grandi & Molendi 2002), it is sufficient for estimating the 
magnitude of the systematic error. Fitting observations with 
eq. (12) gives 7 = 1.46 ± 0.06 for non-cooling flow clusters, 
and 7 = 1.20 ± 0.06 for cooling flow clusters (De Grandi & 
Molendi 2002). With 7 = 1.46 one gets a systematic shifts 
of the peculiar velocity of 20% (for 7 = 1.20 one gets 6%), 
and for the theoretically extreme case of an adiabatic gas 
with 7 = 5/3 the error becomes about 30%. A simple fit to 
the error in percent of the peculiar velocity as a function 
of 7, which is accurate for 7 > 1.2, gives: error (in %) = 
53(7 - 1) - 5. 

Cluster merging induces bulk motion within the cluster 
as large as 500 km/sec (see e.g. Sunyaev et al. (2003)), and 
the corresponding mis-estimate of the peculiar velocity is 
about 10% (Haehnelt & Tegmark 1980), or up to about 100 
km/sec (Holder 2002). 



4 ONION PEELING A SZ CLUSTER 

We have seen above that the SZ determined peculiar velocity 
of realistic galaxy clusters may be systematically wrong by 
10 — 20%, and that X-ray observations cannot be expected 
to solve this problem because they are observationally ex- 
pensive. The natural question is then: what is the solution? 



© 2002 RAS, MNRAS 000, 1-4 



4 Steen H. Hansen 



For X-ray observations the similar problem is solved by 
the technique of deprojection. This heuristically corresponds 
to peeling an onion. The outermost layer is observed and 
analysed. Then the next layer is analysed while subtracting 
the signal from the outer layer. The error-bars will increase 
for the inner layers, but for X-rays this turns out not to 
be disastrous. That is because the contribution from the 
outer layers is significantly smaller than the contribution 
from the inner layers, since the emission drops fast with 
radius e ~ T 1/2 n 2 ~ r Cx with ( x ~ -4.2. 

Similarly, for SZ observations the relevant quantity is 
the local Compton parameter, which goes like n e T e ~ r c ", 
with £ H ~ —2.5. Thus, the error-bars should remain under 
control even when subtracting contribution to the intensity 
from outer layers. Naturally, the fact that 2.5 < 4.2 is the 
reason why SZ observations are more suitable than X-ray 
observations for measuring the outer cluster region. 

An advantage for SZ deprojection is that the peculiar 
velocity is the same for the entire cluster when we can ig- 
nore the local gas motion (see however Holder (2002)). That 
means that while normally one would expect to have 3 SZ 
parameters (y,T e ,v p ) for each of the N radial bins, giv- 
ing 3JV parameters, then we really only have 2N + 1 free 
parameters. This is important because the major parame- 
ter degeneracy for SZ observations is between the tempera- 
ture and peculiar velocity. Thus, the deprojection itself may 
help reduce the error-bars on the derived parameters. It is 
worth mentioning that spatially resolved SZ observations 
will provide both temperature and electron density profiles, 
so combined with accurate X-ray observations one can di- 
rectly measure the radial dependence of dumpiness, which 
is a non-trivial measure of the merger history. 

Such SZ deprojection of the observed intensity clearly 
demands that the future SZ observations should both have 
good spectral coverage and good spatial resolution. That is 
a non-trivial observational challenge. 



5 CONCLUSION 

We have shown that due to the non-trivial temperature 
structure of galaxy clusters, the SZ derived peculiar velo- 
city will be systematically shifted by 10 — 20%. The Comp- 
ton weighted electron temperature is, however, derived ac- 
curately. For future all-sky SZ surveys one cannot rely on 
the observationally expensive X-ray observations to remove 
this systematic error, but one should instead reach for suffi- 
cient angular resolution to perform a deprojection in the SZ 
spectra. 
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